Combined potential and spin impurity scattering in cuprates 
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We present a theory of combined nonmagnetic and magnetic impurity scattering in anisotropic 
superconductors accounting for the momentum-dependent impurity potential. Applying the model 
to the d-wave superconducting state, we obtain a quantitative agreement with the initial suppression 
of the critical temperature due to Zn and Ni substitutions as well as electron irradiation defects in 
the cuprates. We suggest, that the unequal pair-breaking effect of Zn and Ni may be related to a 
different nature of the magnetic moments induced by these impurities. 



Impurities offer a useful experimental probe of the fundamental properties of high temperature superconductors. 
Controlled substitutions of 3d transition metals (Zn, Ni) provide indirect information on the nature of the pairing 
mechanism since they affect the physical properties of both the superconducting and normal state. The nominally 
nonmagnetic Zn (3d 10 , S = 0) as well as magnetic Ni (3d 8 , S = 1) atoms-reside in the magnetically active in-plane Cu 
(3d 9 , S = 1/2) sites. Macroscopic susceptibility and NMR measurementsEHj of Y Ba2Cu^Oj^s (Y-123) superconductor 
indicate the impurity-induced magnetic moments of 0.86/is for Zn, and 1.9fiB for Ni in the underdoped compound, 
which decrease with hole doping to 0.36fiB/Zn and LQ^s/Ni in the optimally doped system. In La,2- x Sr x CuO\ 
(La-214) the same impurity substitutions lead to magnetic properties corresponding to local magnetic moments of 
l.Q^B /Zn and 0.6/is /Nin Magnetism generated by Zn and Ni atoms is, however, different in nature™ JAJhereas, 
the Zn-induced moments reside in the vicinity of the nonmagnetic impurity on the neighboring Cu sitespu'Q the Ni 
substitution yiclds-tlie magnetic moments which partially screen S = 1 impurity spin and result in a more localized 
5=1/2 moment.ElGI In the cuprates, where the strong Coulomb interactions are present, a creation of a magnetic 
momcat-around any lattice defect cannot be excluded. Particularly, a removal of the oxygen atom from the CuOi 
planejaiiil where the Cu (3c?) spins and O (2p) holes are strongly coupled, should give rise to a local uncompensated 
magnetic moment, that is it should develop a magnetically active oxygen vacancy defect. 

A theoretical analysis of the impurity influence on the d-wave superconducting state based on a standard isotropic 
(s-wave) scattering approach indicates a significantly stronger suppression of superconductivity than observed in the 



(s- wave J scattering approach indicates a significantly stronger suppression ot superconductivity than obsepfee m the 
cuprates Jl3 This discrepancy can be reconciled within a simple model of anisotropic impurity scatteringpj'tj which 
takes also thfijipn s-wave scattering channels into account. Evaluation of the critical temperature^ and the upper 
critical fielctJ'Eil dependences on the scattering rate showed their robustness to impurity scattering in the d-wave 
channel and led to a quantitative agreement with the experimental data. A strong support to these early theoretical 
predictions is provided by a recent experimenfO with the electron irradiation created oxygen defects in Y-123 which 
confirms the extended nature of the scattering centers and the ability of the model to account for both the T c and H C2 
suppression by the disorder. Despite its preliminary success this approach neglects the magnetic scattering processes 
which however should be present at least for Zn and Ni impurity substitutions. Obviously, the spin- flip scattering 
leads to an additional destruction of a singlet d-wave superconducting state and, if included in our model, could 
possibly shift its predictions beyond the experimental rang£- lr _rn this paper, we discuss this issue in more detail 
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by generalizing the effective anisotropic impurity scatteringaLj to a combined potential and magnetic scattering. 
We analyze the-jioitial T c sw«pression in Y Ba2Cu^,0'j^s, La 2 - x Sr x CuOi, Bi2Sr 2 CaCu20s+s (Bi-2212) compounds 
induced by-ZptaTHJ and NiEHEj substitution as well as by the in-plane oxygen vacancy defect created by the electron 
irradiation.Ertil Discussing the T c alone cannot give a final test of the impurity scattering model. Its applicability can 
be ultimately verified by a thorough analysis of the thermodynamic and transport properties as has been described 
for isotropic scattering in Refs. 26-28. Such a comprehensive approach is beyond the scope of this brief report in 
which we focus on a possible role of the extended magnetic scattering in the suppression of the critical temperature. 
An elaborated study of the presented model will be given in a separate paper. 

We consider randomly distributed impurities at low concentration interacting with conduction electrons through a 
potential u(k, k') = v(k, k') + J(k, k')S(7, where S is a classical spin representing the impurity and a is the electron 
spin density. The anisotropic superconducting state is determined by the order parameter A(k) = Ae(k), where e (k) 
is a momentum-dependent function. We normalize e(k) in this way that (e 2 ) = 1, where < ... >= J FS dSkn(k)(...) 
denotes the average value over the Fermi surface (FS) momenta, and n(k) is the angle resolved FS density of states 
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normalized to unity, i.e. J FS dSkn(h) = 1. For the (d x 2_ y 2 + s)-wave state e(k) = (cos2(f) + s) J ( (cos2<fi + s) 2 \ in a 
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polar angle notation where the d x 2_ y 2-w&ve state corresponds to s = 0. Taking the electron-impurity scattering within 
the Born approximation and neglecting the impurity-impurity interaction,!^ the diagcm&l-and off-diagonal corrections 
to the Green's function averaged over the impurity positions and spin directions readE3'E3 



Co (k) = to + nn 
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where to — nT(2n + f ) is the Matsubara frequency (T is temperature and n is an integer number), Nq is the single- 
spin density of states at the Fermi level, n$ is the impurity (jdef^t) concentration, and all the wave vectors are 
restricted to the Fermi surface. We generalize our earlier modell3'Ej to include magnetic scattering and assume that 
the momentum-dependent potential terms in Eqs. (H), (0) are separable and given by 



(k,k') 
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J 2 (k,k') = J 2 + J 2 . 9 (k) 5 (k') 
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where vq (t>i), Jo (Ji) are isotropic (anisotropic) scattering amplitudes and / (k), g (k) are the momentum-dependent 
anisotropy functions in the nonmagnetic and magnetic scattering channel, respectively. 

We assume that the overall scattering rate is determined by the isotropic components and to have it well defined 
impose the constraints 
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t2 < t2 



(/) = 0, 
(a) = 0, 



</ 2 > = 1 
(a 2 ) = 1 
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Note, that such a choice of the anisotropy does not change the normal state properties as the diagonal part of the 
self-energy in the limit of A — > is given by to (k) A=0 = uj + imiNo (vq + J 2 S (S + 1)) sgn (to) and depends only on 
the isotropic scattering rates. 

Inclusion of a momentum-dependent model impurity potential in the magnetic channel is motivated by the existence 
of extended magnetic moments associated with the impurity,Enij and the anisotropy of the nonmagnetic scattering 
channel .simulates the anisotropy of the crystal lattice. For Ni impurity which leads to a screened localizecLroagnetic 
moment&Q we expect Jx/Jo -C 1, while for broadly distributed over Cu sites Zn- induced magnetic momentETEJiZI both 
isotropic and anisotropic scattering channels should be comparable in magnitude, i.e. Ji / Jo ~ 1 ■ Contrary to the NMR 
measurementaSa muon spin rotation (/iSR) experimentsEU report no evidence for additional Zn-induced moments in 
underdoped Y-123. Yet, one cannot exclude the possibility that the impurity-generated magnetic moments exist but 
are weakly coupled to the spin system of the Cu0 2 planes. Therefore, the following assumption about the relative 
amount of magnetic scattering seems reasonable: Ji/ui < 1, Jo/vq < 1. __ 

Taking a separable pair potential V (k, k') = — V>e (k) e (k'), Vq > 0, and following a standard procedureEj'El we 
obtain the critical temperature, T c , equation 
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where T Co is the critical temperature in the absence of impurities, ij) (z) is the digamma function, 
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and To = nniNoVQ, Tx = miiNovf, Go — nniNoJoS (S + 1), and Gx — nriiNo J 2 S (S + 1) are the intrinsic scat- 
tering rates. Given strongly coupled spin and charge dynamics in the copper-oxygen planes the impurity potential 
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encountered by the electron in the spin scattering channel should follow the one in the nonmagnetic channel, that is 
g (k) = ±f (k), which simplifies Eq. (0) to 
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Comparing the above equation with the one for nonmagnetic scattering (Eq. (24) of Ref. 13) we notice that there is 
only one additional term in Eq. (||) reflecting an extra T c suppression due to the spin-flip scattering determined tw 
the exchange rate 2Go. The remaining terms of Eq. (||) correspond to the ones obtained for potential scatteringli^l 
with some new effective scattering rates being combinations of the scattering rates in the magnetic and nonmagnetic 
channels: Tq = Tq + Go, Ti = Ti — G\. We note, that the anisotropy of the impurity potential has no effect on the 
T c when T± — 0, i.e. when I\ = G±, which should not be the case of the cuprates, where most probably the exchange 
interaction between the impurity spin and the Cu-0 spin system is small, G\/T\ < 1. The role of anisotropy in the 
impurity potential is concurrently expressed by the dimensionless parameter (ef) (0 < (ef) < 1) representing the 
interplay between the pair potential and the anisotropic part of the impurity potential. Again, for (ef) |-^j= we obtain 
a regular T c suppression due to the isotropic combined nonmagnetic and magnetic impuxity scattering^ which in the 
case of a conventional s-wave superconductor reduces to the well known relationshipBSEEl 

In the quantitative analysis of the experimental data we consider the T c reduction in the limit of low impurity 
concentration n., — > which is given by the initial slope 



dr 4 



>M 



(9) 



where for the sake of brevity we have introduced x = 1 — ( e ) , <j>A = ( e /) Ii/ITo (anisotropy factor), and 4>m = 2Go/r 
(magnetic factor). The isotropic part of the impurity potential Tq determines the residual resistivity at the zero 
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within the ra: 



fiflcm, where u> p i is the in-plane plasma frequency in eV. We fix ui p i 



10.18 x 10- 8^T C0 ^ 

12 l.leV < uipi < lAeV for Y-123, and estimate it from the zero temperature penetration depth 



measurementSc^ as 0.84ey and 0.9eV for La-214 and Bi-221^j«§Bectiyely. There is still no C Qiupl ete experimental 
agreement on the impurity-induced resistivity in the cuprates.tSEjOoEil However, some resultalaO suggest that Zn 
and Ni lead to a comparable effect in the normal state, which would mean that f o (Zn) ~ Tq (Ni). The residual 
resistivity po is used to express the initial T c suppression in the following form convenient for the discussion of the 
experimental data 
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We find that in the d-wave state (x = 1) the initial T c suppression (Eq. (|10[)) is determined only by the dimensionless 
anisotropy factor <pA- The influence of the spin- flip scattering is present through the renormalized scattering rates 
in the isotropic and anisotropic scattering channels: Tq and Ti. However, the magnetic factor <pM has no effect on 
the critical temperature and we may say that the d-wave system [becomes insensitive to magnetic scattering as the 
solutions of Eq. (^|) correspond to those for nonmagnetic scatteringp] with renormalized parameters. In Tabs. I-III we 
give the anisotropy factors cj>A needed to account for the experimental data of Zn, Ni and oxygen vacancy-induced T c 
suppression in Y-123, La-214 and Bi-2212 compounds.L3 Worth noting are generally significantly lower values of 4>a for 
Zn substitution than for Ni. This feature can be attributed to a different nature of the magnetic moments associated 
with these impurities. While the Zn-induced moment is broadly distributed over neighboring Cu sites (G\/Gq ~ 1), 
the one created by Ni is screened on a larger distance and is more localized (G\/Gq <C 1). Therefore, for comparable 
amounts of isotropic scattering by these impurities, Tq = To + Go, the anisotropy factor ipA ~ (T\ — G\) / (Tq + Go) of 
Ni scattering potential, ~ T\/ {Tq + Go), should exceed the one of Zn impurity ~ (Tj — Go) / (ro + Go). Given large 
(f>A values for the oxygen vacancy we may also infer that the magnetic moments formed around this defect are rather 
localized or the anisotropic scattering takes place in the d-wave channel, i.e. (ef) 2 ~ 1. Still, to decide definitely 
about the role of spin scattering in high T c materials with simple defects a quantitative estimate of the ratios of the 
magnetic to nonmagnetic scattering rates Gq/Tq and G\/T\ is required. Of some help here may be a generalization 
of the H C2 critical field analysis to magnetic impurity scattering. So far, an application of the nonmagnetic scattering 
expression for the reduced H C2 slope&3ii3 in the interpretation of the critical field of the electron irradiated Y-123 has 
confirmed the anisotropy parameters of the model determined from T c suppression in this compound.tll 
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The orthorhombic distortion of the crystal lattice in the Y-123 superconductor leads to the (d+s)-wave supercon- 
ductivity in this compoundO The photoemission, tunneling and thermodynamic measurements put an upper bound 
of ~ 10% on the relative size of the s-wave component in the (d+s)-wave superconducting state.c3 If we define a per 
cent fraction A s % of the s-wave component as a ratio of a minimal gap to its maximum value, then the condition 
A s% = 10% is equivalent to x = 0.976 ((e) J = 0.0241). It is a matter of a straightforward calculation to show, that 
cf)A for a given superconducting state is related to the one of the d-wave state <fi A through (pA = 4>\ — (1 — x) (1 ~ 4*m)- 
In what follows, for x = 0.976 we get merely a slight change of the anisotropy factor (f>A = 4>a~~ 0.024 (1 — 4>m) which 
for different amounts of the magnetic scattering rate reads: §a = <\>\ — 0.024 for tpM — (Go = 0); 4>a = 4>\ — 0.008 
for <j>M — 0.666 (Gq/To ==..0.5); <j>A — 4>a f° r $ M ~ •"■ (Go/Fq = 1). Our approach reduces to the isotropic impurity 
scattering approximation^ in the limit of 4>a = 0. Therefore, given the nonzero 4>a values in Tabs. I-III and taking 
the above corrections for the (d+s)-wave state into account, we note that the experimental data cannot be explained 
by applying the isotropic impurity scattering approach to the (d+s)-wave superconducting state with an appropriate 
(for the cuprates) amount of the s-wave component. 

Concluding, we find that even in the presence of the spin- flip processes the scenario of anisotropic impurity scattering 
accounts quantitatively for the experimentally observed T c suppression due to Zn, Ni and oxygen vacancy scattering. 
Moreover, a different partition of the magnetic scattering rates into isotropic and anisotropic scattering channels can 
be used in the interpretation of the weak Ni-induced and relatively strong Zn-induced pair-breaking. Our result also 
shows that the isotropic impurity scattering approximation does not explain quantitatively the initial suppression of 
the critical temperature in the cuprates. Finally, we note that a significant sample-dependence of the discussed data 
with the discrepancy even as high as 50% calls for further more controlled measurements of the impurity-induced T c 
suppression in the high temperature superconductors. 

The authors would like to thank Professor Peter Fulde for helpful comments and hospitality at Max-Planck-Institut 
fur Physik komplexer Systeme in Dresden where the manuscript was prepared. 
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TABLE I. Y-123: 4>a anisotropy factor of the impurity potential reproducing the T c suppression in the Y-123 compound 
within the d-wave superconductivity scenario, i.e. x — 1- (1-leV < lo p i < lAeV; o«.=overdoped, op. =optimally doped, 
un. =underdoped) 



defect 


sample 




anisotropy factor : 4>a 


Zn 
impurity 


(ot>.) single crystal 1 " 
(ov.) single crystal 19 ' 22 
(op.) thin film 23 
(op.) thin film 24 


-0.674 
-0.57 
-0.241 
-0.520 


0.093 < 4>a < 0.440 
0.233 < 4>a < 0.526 
0.676 < 4>a < 0.780 
0.300 < 4> A < 0.568 


Ni 
impurity 


(ov.) ceramic sample 22 
(ov.) film 25 


-0.333 
from (-0.063) to (-0.044) 


0.552 < 4> A < 0.723 
0.915 < 4>a < 0.963 


oxygen 
vacancy 


(ov.) single crystal 9 
(ov.) film s 


-0.30 ±0.04 
-0.187 


0.542 < 4> A < 0.784 
0.748 < 4>a < 0.845 


TABLE II. La-214: <j>A anisotropy factor of the impurity potential reproducing the T c suppression in the La-214 compound 
within the d-wave superconductivity scenario, i.e. x = 1- i^pl = 0.84eF; notation as in Tab. I.) 


defect 


sample 




anisotropy factor : 4>a 


Zn 
impurity 


(ov.) single crystal 19 
(un.) film 20 ' 21 


-0.37 
-0.233 


0.146 
0.462 


oxygen 
vacancy 


(un.) film 10 


-0.127 


0.707 


TABLE III. Bi-2212: cf>A anisotropy factor of the impurity potential reproducing the T c suppression in the Bi-2212 compound 
within the d-wave superconductivity scenario, i.e. x = 1- i^pi = 0.9eV; notation as in Tab. I.) 


defect 


sample 




anisotropy factor : 4>a 


oxygen 
vacancy 


(ov.) single crystal 11 


-0.28 


0.437 
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